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Tables of Energy Deposition Distributions in Aluminum and Copper Irradiated 
f * 
by Point-Monodirectional Electron Beams with Energies from 1 to 60 MeV 


M. J. Berger and S. M. Seltzer 
National Bureau of Standards 
Washington, D.C. 20234 


This report presents tables of energy deposition distributions in 
aluminum and copper media irradiated by monoenergetic point-monodirectional 
electron beams. The distributions are given as functions of the depth in 
the medium and of the radial distance from the incident beam. Results are 
given for 7 beam energies between 1 MeV and 60 MeV. As Shown earlier in 
National Bureau of Standards Report NBSIR 82-2451 (1982), the tabulated 
results from monoenergetic, point-monodirectional sources can be used, by 
Superposition, to obtain spatial distributions of the energy imparted to 
the medium by parallel beams with arbitrary spectra and finite arbitrary 
cross sections. 


Key words: absorbed-dose distribution; aluminum; copper; electron; point- 
monodirectional beam; superposition. . 
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1. INTRODUCTION 


The spatial distribution of energy imparted to a medium irradiated by an electron 
beam depends on the electron stopping power and on the multiple scattering of the 
electrons in the medium. Even when the incident beam is narrowly collimated (point- 
monodirectional), the distribution of imparted energy spreads laterally with depth in 
the medium, as the result of multiple scattering. The higher the atomic number of 
the medium, the greater are the effects of multiple Scattering. 

In a previous report ciy! we have presented comprehensive tables of the spatial 
distribution of imparted energy in a water medium irradiated by point-monodirectional 
monoenergetic electron beams with energies from 1 to 60 MeV. It was also shown in 
reference [1] how these "elementary" distributions can be combined by superposition 
to obtain the spatial distribution of imparted energy from parallel electron beams 
with finite cross sections of arbitrary shape. Of course one can also use weighted 
sums of elementary distributions for monoenergetic sources to obtain results for 
incident electron beams with arbitrary spectra. 

The results for water were generated mainly with a view towards applications in 
biomedical dosimetry, for example in radiation treatment planning. In the present 
report the calculations are extended to other media of higher atomic number (aluminum 
and copper) to provide information of interest for radiation damage studies in terms 
of imparted energy. 


2. CONDITIONS OF CALCULATION 


The calculations were, as in reference [1], carried out with the Monte Carlo 
transport program ETRAN [2]. A point-monodirectional electron beam with kinetic 
energy Ue, was assumed to be incident along the z-axis onto a semi-infinite medium 
that occupies the region z2O. The electron tracks simulated by random sampling 
were analyzed to obtain an estimate of the elementary distribution function F(p,z). 
This distribution is defined such that F(p,z)dodz represents the average fraction 
of the incident beam energy that is deposited at depths between z and z+ dz and 
at radial distances between op and op+dpo from the z-axis. This distribution is 


normalized such that 
f az f F(gezidou=ile ans (1) 
) (0) 


The Monte Carlo results were actually obtained, and were tabulated, in the form 
of the following histogram distributions: 


‘Numbers in brackets indicate literature references at the end of this paper. 
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a) Distribution of imparted energy in depth intervals, 
Zz, 90 
H. = 100 f dz fae" Fil(ou Ze) mee (2) 
n 
Zn-1 (0) 


b) Cumulative radial distribution of imparted energy in indicated depth 
intervals, 


An 


— 


p 
dz f do" F(p' 52) 
] (0) 


eae 
¥, (0) = 100 riz wes, TD, ATO KS : (3) 
f dz f do'F(p',z) 
Zy-1 (0) 
The boundaries of the depth intervals were chosen such that 
200% nr/20, NP =~Osads Carre (4) 


where in is the c.s.d.a. range (rectified pathlength computed in the continuous- 
slowing-down approximation) at the incident-beam energy To: Similarly, tabular 
values of V, (0) are given for p-values that are specified fractions of ro: 
distributions of imparted energy are thereby "scaled" in the sense that their 
dependence on the incident beam energy Ue is minimized, which facilitates interpola- 
tion with respect to T5: 


The 


3. ELECTRON RANGES 


The Monte Carlo calculations for this report were done some time ago, with the 
best cross sections then available for the electron energy loss per unit pathlength. 
Since then, the collision stopping power (mean energy loss per unit pathlength due to 
excitations and ionizations) and the radiative stopping power (mean energy loss per 
unit pathlength due to the emission of bremsstrahlung) have been revised [3]. For the 
collision stopping power, more up-to-date experimental values were used for the mean 
excitation energy (the key parameter of the Bethe stopping power formula). Moreover, 
the density-effect correction (reduction of energy loss due to the polarization of the 
medium) was improved. The universal approximation formula of Sternheimer and Peierls 
[4] was no longer used. Instead, the density-effect was evaluated according to the 
original method of Sternheimer [5], using up-to-date values of the mean excitation 
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energy and of the atomic binding energies needed for the modeling of the dispersion 
properties of the medium. In regard to the radiative stopping power, the principal 
changes consist of the use of the new theoretical cross sections of Pratt et al. [6] 
below 2 MeV, the re-evaluation of the cross sections of Davies, Bethe, Maximon, and 
Olsen [7] above 50 MeV, and the use of interpolation with respect to electron energy 
to cover the region from 2 to 50 MeV. 

Fortunately the use of improved cross sections resulted in relatively small 
changes in the stopping powers for aluminum and copper compared to those used in the 
Monte Carlo calculations, and the corresponding changes of the c.s.d.a. ranges, shown 
in table 1, are equally small. It is expected to be a good approximation to assume 
that the scaled spatial distributions of imparted energy, expressed in terms of the 
variables p/ry and z/Y5» are not affected by the change of cross sections. The 
appropriate adjustment to the new cross sections can then be made simply by using the 
new range values when converting the spatial distributions from the scaled form into 
unscaled distributions as functions of p and z. 


4. RESULTS 


Spatial distributions of imparted energy from point monodirectional electron 
beams, are given for 7 beam energies from 1 MeV to 60 MeV. The results at each 
energy are based on a sample of 10,000 tracks of primary electrons as well as the 
tracks of all generations of secondary electrons (including those set in motion by. 
bremsstrahlung photons). Tables 2A and 2B give the histogram distributions Hy 
with respect to depth, for aluminum and copper. Tables 3A and 3B give the cumulative 
histogram distributions V, (0) with respect to radial distance from the incident 
beam, in various depth intervals, for aluminum and copper. 
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Table 1. Electron ranges in aluminum and copper 


ES |e ete 


Percent Percent 
Difference Difference 


20.29 

16.24 

10.43 
6.183 
3.408 
1.384 
0.6367 


* 
Ranges are computed in continuous-slowing-down approximation and represent 
rectified pathlengths traveled by electrons in the course of slowing down to 
rest. 


a 
Old: used in Monte Carlo calculations. 


PNew: From reference [3]. 


TABLE 2A. Distribution of energy deposition in aluminum as a function 
of depth. The numbers in the column headed "Scaled depth 
interval" are the depth boundaries in units of the css.d.a. 
range ro. The other columns give the quantity H, as 
defined by Eq (2), z<.e., the percentage of the beam 
energy that is deposited in the indicated depth interval. 


Energy Deposition as a Function of Depth in Aluminum 


Beam energy, MeV 


Scaled depth 60 40 20 10 5 2 1 
interval 
0.0 - 0.05 2.97 3.38 3.91 4.29 4.61 4.91 5.04 
0.05 - 0.10 3.16 3.58 4.18 4.68 5.20 5.97 6.45 
0.10 - 0.15 3.23 3.73 4.35 5.03 5.88 7.06 7.81 
0.15 - 0.20 3.35 3.76 4.51 5.21 6.43 8.15 8.74 
0.20 - 0.25 3.47 3.88 4.58 5.61 6.95 8.78 9.17 
0.25 - 0.30 3.63 3.98 4.76 6.01 7.64 9.00 9.26 
0.30 - 0.35 3.64 4.07 5.05 6.35 7.99 8.89 8.88 
0.35 - 0.40 3.69 4.14 5.09 6.65 8.06 8.41 8.31 
0.40 - 0.45 3.76 4.16 5.33 6.97 7.76 7.60 7.60 
0.45 - 0.50 3.75 4.19 5.39 6.89 7.31 6.82 6.57 
0.50 - 0.55 3.75 4.27 5.48 6.70 6.58 5.92 5.42 
0.55 - 0.60 3.71 4.25 5.49 6.17 5.81 4.77 4.27 
0.60 - 0.65 3.64 4.20 5.31 5.64 4.92 3.62 3.15 
0.65 - 0.70 3.54 4.07 5.01 4.91 3.89 2.64 2.07 
0.70 - 0.75 3.45 3.82 4.52 3.95 2.80 1.78 1.31 
0.75 - 0.80 3.25 3.60 3.95 3.02 1.87 1.04 0.75 
0.80 - 0.85 3.03 3.38 3.25 2.12 1.06 0.53 0.40 
0.85 - 0.90 2.84 2.99 2.55 1.35 0.54 0.21 0.14 
0.390 - 0.95 2.55 2.56 1.82 0.71 0.24 0.06 0.04 
0.95 - 1.00 2.26 2.11 1.21 0.30 0.07 0.02 0.01 
1.00 - 1.05 1.98 1.69 0.71 0.13 0.03 0.00 0.00 
1.05 - 1.10 1.69 1.33 0.40 0.06 0.02 0.01 0.00 
1.10 - 1.15 1.43 0.98 0.23 0.05 0.01 0.00 0.00 
1.15 - 1.20 1.16 0.72 0.16 0.05 0.02 0.00 0.00 
1.20 - 1.25 0.99 0.52 0.14 0.05 0.01 0.01 0.00 
1.25 - 1.30 0.83 0.40 0.13 0.04 0.01 0.00 0.00 
50 e= Hes 0.72 0.33 0.13 0.03 0.01 0.00 0.00 
1.35 - 1.40 0.60 0.32 0.12 0.03 0.01 0.00 0.00 
1.40 - 1.45 0.56 0.31 0.10 0.03 0.01 0.00 0.00 
1.45 - 1.50 0.53 0.29 0.10 0.02 0.01 0.00 0.0 
1.50 - 1.55 0.51 0.28 0.10 0.02 0.01 0.00 0.00 
1.55 - 1.60 0.48 0.28 0.09 0.02 0.01 0.00 0.00 
1.60 - 1.65 0.48 0.29 0.07 0.02 0.00 0.00 0.00 
1.65 - 1.70 0.47 0.26 0.08 0.03 0.01 0.00 0.00 
1.70 - 1.75 0.44 0.27 0.07 0.02 0.00 0.00 0.0 
1.75 - 1.80 0.42 0.28 0.07 0.02 0.01 0.00 0.0 
1.80 - 1.85 0.39 0.26 0.07 0.02 0.00 0.00 0.00 
1.85 - 1.90 0.39 0.24 0.07 0.02 0.00 0.00 0.00 
1.90 - 1.95 0.39 0.23 0.07 0.02 0.00 0.00 0.0 
1.95 - 2.00 0.34 0.23 0.06 0.02 0.00 0.00 0.0 


TABLE 2B. Distribution of energy deposition in copper as a function 
of depth. The numbers in the column headed "Scaled depth 
interval" are the depth boundaries in units of the c.s.d.a. 
range Yo. The other columns give the quantity Hy, as 
defined by Eq (2), z.e., the percentage of the beam energy 
that is deposited in the indicated depth interval. 


Energy Deposition as a Function of Depth in Copper 


Beam energy, MeV 


Scaled depth 60 40 20 10 5 2 1 
interval 
0.0 - 0.05 2.235 2.66 3.48 4.41 5.62 7.12 7.61 
0.05 - 0.10 2.47 2.95 3.90 5.31 7.37 9.70 10.73 
0.10 - 0.15 2.64 3.13 4.32 6.16 8.68 11.33 12.01 
0.15 - 0.20 2.84 3.42 4.65 6.84 9.56 11.43 11.63 
0.20 - 0.25 3.02 3.60 5.06 7.50 9.67 10.56 10.54 
0.25 - 0.30 3.20 3.75 5.40 7.81 9.32 9.38 9.04 
0.30 - 0.35 3.36 3.96 5.64 7.74 8.55 7.82 7.37 
0.35 - 0.40 3.44 4.04 5.78 7.42 7.62 6.28 5.43 
0.40 - 0.45 3.60 4.16 5.84 6.80 6.38 4.87 4.18 
0.45 - 0.50 3.55 4.13 5.57 6.05 5.16 3.45 2.87 
0.50 - 0.55 3.55 4.11 5.08 5.19 3.94 2.44 1.87 
0.55 - 0.60 3.41 3.83 4.69 4.13 2:85 1.52 1.07 
0.60 - 0.65 3.33 3.60 4.08 3.32 1.88 0.91 0.58 
0.65 - 0.70 3.14 3.42 3.45 2.42 1.18 0.41 0.27 
0.70 - 0.75 2.94 3.08 2.86 1.73 0.65 0.17 0.08 
0.75 - 0.80 2.76 2.78 2.24 1.12 0.33 0.06 0.03 
0.80 - 0.85 2.58 2.41 1.73 0.70 0.14 0.02 0.02 
0.85 - 0.90 2.33 2.01 1.24 0.36 0.08 0.02 0.00 
0.90 - 0.95 2.15 1.72 0.92 0.19 0.03 0.01 0.00 
0.95 - 1.00 1.86 1.44 0.65 0.14 0.04 0.01 0.00 
1.00 - 1.05 1.70 1.25 0.41 0.10 0.03 0.01 0.00 
1.05 - 1.10 1.52 1.06 0.33 0.09 0.03 0.00 0.00 
1.10 - 1.15 1.37 0.93 0.27 £40.08 0.02 0.01 0.00 
1.15 = 1.20 1.19 0.84 0.23 0.10 0.03 0.00 0.00 
1.20 - 1.25 1.03 0.74 0.20 0.08 0.02 0.00 0.00 
1.25 - 1.30 0.94 0.63 0.20 0.07: 0.03 0.00 0.00 
1.30 - 1.35 0.89 0.57 0.19 0.06 0.01 0.00 0.00 
1.35 - 1.40 0.79 0.55 0.17 0.05 0.02 0.00 0.00 
1.40 - 1.45 0.76 0.55 0.17 0.05 0.02 0.00 0.00 
1.45 - 1.50 0.70 0.53 0.19 0.06 0.02 0.00 0.00 
PoBO 2 Wo 0.66 0.52 0.17 0.06 0.02 0.00 0.00 
1.55 - 1.60 0.61 0.47 0.16 0.05 0.01 0.00' 0.00 
1.60 - 1.65 0.60 0.45 0.16 0.05 0.01 0.90 0..00 
1.65 - 1.70 0.60 0.43 0.14 0.05 0.01 0.00 0.0 
1.70 - 1.75 0.59 0.40 0.14 0.04 0.01 0.00 0.00 
1.75 - 1.80 0.55 0.38 0.14 0.03 0.01 0.00 0.00 
1.80 - 1.85 0.53 0.38 0.14 0.04 0.01 0.00 0.00 
1.85 - 1.90 0.55 0.34 0.15 0.04 0.01 0.00 0.00 
1.90 - 1.95 0.52 0.33 0.15 0.03 0.01 0.00 0.00 
1.95 - 2.00 0.48 0.30 0.13 0.02 0.01 0.00 0.00 
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